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ABSTRACT 



We employ raPI effective action techniques to study M = 4 super Yang-Mills, and write 
down the 2PI effective action of the theory. We also supply the evolution equations of 
two-point correlators within the theory. 



1 Introduction 



Our current understanding of phenomena in thermal equilibrium is extensive, but many 
of the processes which give deeper insights into our fundamental understanding of nature 
begin far away from equilibrium. There is abundant experimental data concerning the 
early stages of heavy-ion collisions, which requires the development of a nonequilibrium 
theoretical framework to allow for correct interpretation of the data. Furthermore, we might 
be able to make contact with the study of black hole formation and evaporation (see, for 
example [1]) by using the AdS/CFT correspondence [2]. This widely studied correspondence 
postulates an exact equivalence between string theory on the AdS§ x S 5 background and 
(3 + l)-dimensional Af = 4 super Yang-Mills theory (SYM) (for comprehensive reviews, see 
[3] and [4]). It is the best understood example of a gauge/gravity theory duality, i.e. of 
holography, a groundbreaking hypothesis which says that any gravitational theory should 
have a description in terms of a QFT with no gravity in one less dimension. It is precisely 
this N = 4 SYM which we wish to study further in a nonequilibrium setting. 

The problem of studying nonequilibrium phenomena is two-fold, because we not only 
need to take into account quantum fluctuations, but we also need to deal with a very large 
number of degrees of freedom. Classical statistical field theory simply is not good enough, 
and standard perturbative approaches based on small deviations from equilibrium are not 
applicable: secular, time-dependent terms may appear which invalidate the perturbative 
expansion. For example, in [5] it was argued that for SU{N) Yang-Mills on a sphere, the 
high temperature phase of the theory is intrinsically non-perturbative. In recent years, 
so-called nPI effective action techniques have been developed [6, 7, 8, 9, 10], allowing us 
to use nonperturbative approximations to get a handle on nonequilibrium dynamics, in the 
hope of ensuring non-secular and universal behaviour (meaning that the initial conditions 
do not affect the late-time behaviour). Of particular interest is the precise time evolution 
of quantum fields whose initial state is far from equilibrium. Relevant references relating to 
far-from-equilibrium quantum fields and thermalization include [11, 12, 13, 14, 15, 16, 17, 
18, 19, 20]. For a comprehensive review on progress in nonequilibrium QFT, see [21] and 
the references therein, and the book [22]. 

This paper is arranged as follows. In section 2 we present our results: we give the 
2PI effective action of Af = 4 SYM, and also write down the evolution equations of the 
2-point correlators for the scalars, gluons, fermions and ghosts within the theory, in the 
nonequilibrium realm. 
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2 The nPI Effective Action Approach 



In order to study nonequilibrium quantum field theory, one first needs to specify an initial 
state at some time to, which is usually done by specifying a density matrix po(to) which is 
not a thermal equilibrium density matrix. Equivalently one can specify all initial n-point 
correlation functions, although in practice one often supplies only the lowest correlation 
functions at to- The time evolution of this initial state (i.e. these initial correlation func- 
tions) is then encoded in the functional path integral with classical action S. For example, 
in the case of a real scalar field ip, the nonequilibrium generating functional for correlation 
functions is 

Z[Ji, J 2 , ■ ■ ■ ; P d] = Tr | PD (0)T c e<^ U M*,vM*)*to+ -) J ? {2A) 

where $(x) denote Heisenberg field operators. Tc denotes time-ordering along the time 
path C and in what follows, f = f c dx° f d d x. It turns out that the extension to the 
nonequilibrium realm is done precisely via the introduction of this finite, closed real-time 
contour C, known as the Schwinger-Keldysh contour, given in Figure 1. We call the top part 
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Figure 1: Schwinger-Keldysh contour. 

of the contour, the forward piece, C + , and the bottom backward piece C~ . Time ordering 
along the contour is largely intuitive: we want any time on C~ to be later than any time on 
C + , so we use normal time ordering on the forward piece, and antitemporal time ordering 
along the backward piece. Time integration along the contour is given by 

dx°= / dx° + ' dx°= dx° - / dx°. (2.2) 
Jc Jo,c+ Jt Jo,c+ Jo,c- 

In (2.1) we have included n source terms J\(x), J2(x,y), ■ ■ ■ . Now, the standard pre- 
scription for writing down a 1PI (one particle irreducible) effective action involves intro- 
ducing only the J\ source term when setting up the generating functional, but in order 
to extend this to the so-called "nPI" (n particle irreducible) effective action, we need to 
introduce n source terms. In practice, an equivalence hierarchy exists between nPI effective 
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actions, and means that one can avoid having to calculate the effective action for arbitrarily- 
large n. For a (/-loop approximation, this hierarchy states that all nPI descriptions with 
n < q are equivalent, so only the gPI effective action is required. In our analysis of M = 4 
SYM we will be interested in a two-loop approximation, which amounts to calculating the 
2PI effective action of the theory. In (2.1) this would amount to including only the sources 
J\(x) and J 2 (x,y), so that we can rewrite the two-source generating functional as 

Z[J 1 ,J 2 ; P d] = [dtp+dtp- (<P + \ PD (0)\<P-) r pV^ 5M+/ * Jl( * M * )+l ^^ 

J Jip+ 

(2.3) 

where (p^ are eigenstates of the Heisenberg field operators at initial time, namely <3?(t = 
0, x)]^) = ip ± {x)\(f ± ). It follows from (2.3) that the structure of the nonequilibrium par- 
tition function is very similar to the zero temperature and thermal case, apart from the 
additional piece coming from the initial conditions, and the time ordering along C replacing 
the original time ordering along only C + . In principle, barring the initial conditions, we 
should be able to do our calculations at zero temperature, and easily transform our re- 
sults to fit the nonequilibrium case by introducing the Schwinger-Keldysh contour. Indeed, 
notice that in (2.3) the second integral is basically just the vacuum generating functional 
for connected Green's functions for a scalar field theory with classical action S[<p], in the 
presence of two source terms, 

Z[Ji, J 2 ] = e lW ^ h ^ = Jv<p exp (i S[<p] + jf Ji(xMx) + i jf J 2 (x,y)ip(x)<p(y) 

(2.4) 

Now, in order to obtain the equations of motion of the correlation functions of such a 
simple scalar field theory, thereby fully describing it, one first extracts the effective action F 
from S by performing suitable Legendre transforms. In addition to the Legendre transform 
of the generating functional necessary to obtain the 1PI effective action, we simultaneously 
perform a second Legendre transform to get the required 2PI effective action, namely 

r[4>,G\ = W[Ji,J 2 ] - . . Ji(x) - / — v Mx,y) 

Jx SJi(x) J xy 5J 2 {x,y) 

= W^h]- U{x)J l (x)- \ I J 2 {x,y)<p{x)<P{y)-\TrGJ 2 . (2.5) 

Jx Jxy 

In (2.5), <j>{x) is the macroscopic field (i.e. the field expectation value of tp) and G(x,y) 
is the connected two-point function. An nPI effective action would require n simultaneous 
Legendre transforms of this type. Finally, first order functional derivatives of this effective 
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action with respect to the appropriate correlation functions (in the absence of sources) then 
provide the corresponding equations of motion via the so-called "stationarity conditions" , 

S^.O, Q^J=0. (2.6) 
o<p oG 

For more on the 2PI effective action of scalar and fermion fields, see [6, 7, 8, 21] and the 
references therein. 

2.1 2PI Effective Action of Af = 4 SYM 

N = 4 super Yang-Mills (SYM) has an SU(N) colour gauge symmetry and corresponding 
gauge field A^, and also contains four spinors Aj, where i = 1, ... ,4 transforming under 
the global SU(4) symmetry, and six scalars M m , where m = 1, . . . , 6, transforming under 
SO(6). In order to quantize this theory properly we need to gauge-fix, which we do using the 
Faddeev-Popov procedure. Together with the gauge-fixing term — ^Tr(d fJi A^) 2 and ghost 
term Tr(—fjd^(V^r])) (with ghost fields labelled by i] and anti ghosts by 77), the Af = 4 
SYM action is given by 

SSYM = S° SY M + S SYM, ( 2 - 7 ) 



where 



and 



S°sym = Jjr (^-A,d 2 A» + l -M m d 2 M m + f\ M a^dX ~ ^ 2 r?) , (2.8) 



STym = J Tr^-ig^A^A" - d^A" A») + \g 2 {A^A^ A" - A^A" A») 
- ig {8 ll AI m A^M m - d il M m M m A^) + g 2 {A^M m A»M m - M m A^M m ) 
+ X -g 2 {M m M n M m M n - M m M n M n M m ) - g (xfa^-A^ - Xfa^A^ 

+ \i9 (K\ aj (v m 1 y j M m - Af (a- 1 )^M m A Qj - ~X\X^ (a m \ M m + A^<r m )^M m A^) 
+ ig (d^fjA^ - dPfrnA,,)) . 

(2.9) 

In the above, V M is the covariant derivative, with V^r/ = d^rj + iglA^, 77], and we work in the 
Feynman 'tHooft gauge where £ = 1. The free gluon, scalar, fermion and ghost propagators 
(denoted Dq,Sq,Aq, and Go respectively) are given by 

i 1 

D (x,y) = 4?r 2 ^ _ y y = S o( x ,v) = -G (x,y), \^{x-y) = ia^.d*D (x,y). (2.10) 
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The 2PI effective action of M = 4 SYM is given by 

T[M,A,X,fj,S,D,A,G] = Ssym[M, A, A, fj] 

+ ^TrlnS" 1 + ^TrS^S + ^TrlnD^ 1 + l -TrD^ l D 

- iTrlnA- 1 - iTrA^A - iTrlnG^ 1 - iTrG^G 

+ T 2 [M,A,X,fj,S,D,A,G], (2.11) 

where the S, D, A and G are full scalar, gluon, fermion and ghost propagators respectively, 
and M, A, A, fj are the respective macroscopic fields. The first nine terms above represent 
the 2PI effective action to one loop order, while ^[M, A, A, rj, S, D, A, G] is the higher loop 
contribution. ^[M, A, A, fj, S, D, A, G] is obtained by shifting each of the fields in Ssym 
by the respective macroscopic field (the field expectation value), and using the vertices 
obtained from this shifted action to build the higher loop 2PI diagrams. 

At two-loop order, the fourteen diagrams given in Figure 2 contribute to T 2 in (2.11). 
The wavy lines correspond to gluons, straight lines to scalars, arrowed lines to fermions and 




Figure 2: 2PI 2-loop diagrams contributing to T2[M, A, A, fj, S, D, A, G]. Starting at the top row and going 
from left to right, we label them as S 2 , DS, D 2 , D 3 , S 2 D, AAD, A 2 S, and GGD. 

dashed arrowed lines to ghosts. Vertices with black dots correspond to cubic vertices with 
the shift fields. We will keep only the diagrams at leading order in g (i.e. g 2 ). Then all 
diagrams which depend on the macroscopic fields do not contribute, since they are all either 
0(g 3 ) or 0{g 4 ). Hence, we can exclude all diagrams in Figure 2 which contain black dots. 
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We will label the eight leading order diagrams according to the propagators they contain. 
Thus, starting from the top row in Figure 2 and going from left to right, we have S 2 , DS, 
D 2 , D 3 , S 2 D, AAD, A 2 S, and GGD. Thus, 

t 2 [s, d, a, g]sym = r 5 2 + r DS + r D 2 + r D 3 + r S 2 D + r AAD + r A 2 S + r G g D , 

and 

T 5 2 = -15g 2 (N 3 -N) f S 2 (x,x), 

J X 

Tds = -6g 2 (N 3 -N) [ D£(x,x)S(x,x), 

J X 

T D2 = ± g 2 (N 3 -N) J (Dl(x,x)DS(x,x)-DZ(x,x)DZ(x,x)), 

T D3 = -ig 2 (N 3 -N)l (d;dyD UK (x,y)(D^(x,y)D^(x,y)-D^(x,y)D™(x,y)) 

Jxy 

+8%Dl{x,y) {d}DZ(x,y)ir K (x,y) - ^(x^D^y)) 
+d?D5(x,y) {d»D v K {x,y)D^{x,y) - d»D£(x,y)D''f>(x,y))) , 
T S2D = -6ig 2 (N 3 -N)[ D^(x,y)(d;S(x,y)dyS(x,y)-d;dyS(x,y)S(x,y)), 

Jxy 

Taad = -±ig\N 3 -N)a» a»J A^(x, y)A^(y, x)D^(x, y), 

Jxy 

T A 2 S = -2Aig 2 (N 3 -N)e aP e^f A^(x,y)A^(x,y)S(x,y), 

Jxy 

Tggd = ig\N 3 -N)l d^G(y,x)d^G(x,y)D^(x,y). (2.12) 

Jxy 

We can evaluate the 2PI effective action of M = 4 SYM to 0{g 2 ) by using the free prop- 
agators when evaluating each of the diagrams. Due to the conformal and supersymmetric 
nature of M = 4 SYM, we expect the effective action to vanish, and it does. This is a first 
check that our effective action is indeed correct. 

2.2 Evolution Equations of jV = 4 SYM 

Having determined which 2PI diagrams contribute to the 2PI effective action of N = 4 
SYM, we can use the stationarity conditions to write down the equations of motion for 
each of the fields in our theory [15, 21]. This involves writing down the self energy T,(x,y) 
for each propagator in the usual way (by cutting that propagator line on each of the 2PI 
diagrams at our disposal). As mentioned in the beginning of section 2, barring the initial 
conditions, we should be able to do our calculations in the vacuum, and easily transform 
our results to fit the nonequilibrium case by introducing the Schwinger-Keldysh contour. 



6 



For scalars, the equations of motion are 

[d 2 x -^ s \x;S))F^(x,y) = £ dz° J d 3 z^ p s \x, z)F& (z,y) 

- J* dz°j d 3 zrt;\x,z)p( s \z,y), 

(dl-rt°X s \x;S)) p( s \x,y) = j* dz° j d 3 zZ ( p s \x,z)p( s \z,y). (2.13) 

The superscript (s) refers to scalars (and similarly, throughout this paper, the superscripts 
(gZ)=gluon, (/)=fermion and (g/i)=ghost). 

The integro-differential evolution equations in (2.13) are for the scalar statistical propa- 
gator F^ s \x, y) and the scalar spectral function p^(x, y), which are obtained via a splitting 
of the scalar two-point function [14, 15] 

S c (x,y) = F^ s \x,y) - ^ s \x,y)sign c (x° - y°). (2.14) 

The preference in using F(x, y) and p(x, y) is that they are both real, which makes their 
evolution equations intrinsically more manageable, and more importantly they have handy 
physical interpretations. The spectral function involves the spectrum of the theory while 
the statistical propagator deals with occupation numbers. 

The quantities T,^,\x,y) and T, p s \x,y) arise from a similar splitting of the scalar self- 
energy, namely 

^\x,y) = 5c(x,y)^ s \x) + ^\x,y) - l -X p s \x, y)sign c (x° - y°). (2.15) 

The local term ^^^(x) is part of a generalized "mass" term. (When considering the 
evolution equations for the other fields in M = 4 SYM, we perform a similar splitting of 
the two-point functions and corresponding self-energies, and label them by the superscripts 
(<?/), (/) and (gh)). 

Of course, the scalar evolution equations in (2.13) are quite general. In order for them 
to apply specifically to M = 4 SYM, we need to write down the quantities and 
in the context of this theory. These quantities are related to the corrections to the scalar 
propagator, which we obtain by cutting a scalar line in the two- loop diagrams of Figure 
2. The diagrams in Figure 3 contribute. Being careful to take the contour C into account, 
suppressing adjoint indices and working to 0(g 2 N), we obtain 
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Figure 3: Corrections to the scalar propagator. 



4 s) (x,z) 



-g 2 N \8d^F^(x,z)F^ u (x,z) - 2d^p^(x,z)p^ fXl/ (x, 
+ 4d*F( s \x, z)d*F^ u (x, z) - d*p( s \x, z)d z uP ^ u {x 
+ AdlF {s \x, ^fW""^ z) - d z p {s \x, z)d*p W ^{x 
+ 2F ( - s \x,z)d z d*F^> MU (x,z) - ^p( s \x,z)d^p^ u (x,z) 
- e afj£ fa (4F^ afi (x, z)F^(x, z) - p^ a& (x, z)p^(x, z) 



x,z) 
z) 



+ 4F^ ad (z,x)F^{z,x)-p^ aA (z,x)p^(z,x)) 



(2.16) 



-g 2 N \8 (d*d*F (s) (x,z)p (9l ^ u (x,z) - 2d*d*p( s \x, z)F^ v (x, z) 
+ 4 (dlF i - s \x,z)d z uP W ^{x,z) - d^ s \x,z)d z F^ v (x,z) 

+ (^F( 9 \x,z)d£pW^(x,z)-dZp( 9 \x,z)d£FWi»'(x,z) > ) 
+ 2F^(x,z)d^p^ u (x,z) + 2p^{x,z)d^F^ u (x,z) 

~ 4£ aP £ $a {F U)a& {x, Z)p^(x, Z) + p^ aA ( X , z)F^(x, Z 



- FW a *(z,x)pVW(z,x) - p (f)a *{z,x)F (f W{z,x)^J 



(2.17) 



With analogous definitions for the gluons, fermions and ghosts, we can also write down 
the evolution equations for the statistical propagators and spectral functions of these fields. 
For gluons, the equations are 

^^-StW^s;/))) = jT dz° J d z zY,fl K (x,z)F^ l >^(z,y) 

- J"" dz° Jd*z^f:(x,z)p^(z,y), 

[g^-^ l i K (x-D)) p^(x,y) = J* dz° J d 3 z^(x, z)p^ (z,y). 

(2.18) 
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The fermion equations are 

(ia^ + E^ f \x;A))F^(x,y) = J* dz° J d 3 zE^(x, z)F^ (z,y) 

- j" dz°j d 3 z^. / ,(x,z)p^(z,y), 
^dl + ^ f \x-A)) pW\x,y) = jT dz° J ' d 3 z^. p {x,z)p^\z,y). 



(2.19) 



Finally, the ghost equations are 

^-S( »W( S ;G))fW(i,j/) = £ dz° J d 3 zE^ h \x,z)p^ h \z,y) 

- £ dz° J ' d 3 z^ h \x,z)F^ 9h \z,y), 

(d 2 x -E( 'W( I; G)) P {gh) (x,y) = - £ dz° j ' d 3 zrt p 9h \x,z)p( 9h \z,y). 



(2.20) 



The corresponding Hp 1, ' and Y,^ 1 ^^^ 9 ^ for each of the equations above are given 
in the Appendix. 



3 Discussion and conclusions 

In this paper we used the nPI effective action approach to write down the 2PI effective 
action of AT = 4 SYM. We then wrote down the evolution equations for the two-point 
correlators of scalars, gluons, fermions and ghosts in the theory. A particularly pleasing 
property of N = 4 SYM is that it is a finite theory, so no renormalization is required. nPI 
effective actions enable us to set up a very efficient nonperturbative approximation scheme 
for nonequilibrium QFT, in the hope that we may bypass the usual problems of secularity 
and non-universality experienced by the standard perturbative approaches. It would be 
interesting to see the dependence of our results on gauge fixing, since the 2PI effective 
action in general has a gauge dependence [9] . 

An important step in understanding the nonequilibrium dynamics of quantum fields 
is to understand how systems which are initially far from equilibrium approach thermal 
equilibrium at late times. We want to understand thermalization in QFTs which have a 
holographic dual, in particular M = 4 SYM. Thus, a first step in extending our result 
would be to consider various initial conditions and solve the equations of motion, thereby 
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allowing one to explore thermalization in this theory. Arguably the simplest possibility is to 
consider gaussian initial conditions and solve the equations of motion numerically (any such 
solution would have to be numerical, due to the extremely nonlinear and coupled nature of 
the evolution equations). One could them move on to more complicated initial conditions, 
since physical initial conditions may not be gaussian. In particular, since N = 4 SYM is 
a super symmetric conformal field theory, supersymmetric initial conditions may simplify 
things. Exploring thermalization in this theory and comparing it to that in QCD would 
potentially give insights as to why the RHIC data seems to be well-described by the strongly 
coupled M = 4 theory. It would also be interesting to investigate the implications for black 
hole formation and evaporation, since, through holography, the process of thermalization 
is expected to be mapped to horizon formation on the gravitational side. In the context 
of holography, a similar formulation was developed in [23, 24, 25], but instead for gauge 
invariant operators without the 2PI technique. 
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A Appendix 

The corrections to the gluon propagator in N = 4 SYM are given by the six diagrams in 
Figure 4. Thus, we have 
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N cr p .al- \8F^ ap (x, z)F ( -^(z, x) + 2p^ ap \x, z)p^ Kl3 {z, x) 
+12d p F^ (x, z)df>F^ (x, z) - 3d p p^ (x, £)<9f p (s) (x, z) 
-UdZdfF^ (x, z)FW (x, z) + Zd p x dlp [s) (x, z)pW (x, z) 
-2diF^ h \x,z)d p F^ h \z,x) - ^dlp ( - 9h \ Xl z)d p p i9h \z,x) 

2^F(^(i,^( j| W(i,z) - ld*d z ppW p t(x,z)p( 9l ^(x,z 



+4 
+4 
-2 
+4 
+4 
-2 
-2 
-2 
+8 
-4 
-4 
+2 
-4 
+2 
-4 
+2 
+2 



2 

2d*d z pFW p P(x,z)FWrt{x,z) - ^d*d z ppW p P(x,z)p( 9l ^(x,z) 
2^/M^(i,z)FW rf (i,z) - ^d%d z pp( 9l ^(x,z)p < > 9l W(x,z) 
2d*d z pF( 9l W(x,z)FW p t{x,z) - ld*d z pp Wp(3 (x,z)p {9l)p Z(x,z) 



2d p dlF^ pK (x,z)F^(x,z 



2d p dlF^{x,z)Flf{x,z 
2d p d^F^ l \x,z)F^ 1 ^ [X, z 
2d p dlF( 9 J\x,z)F( 9l ^(x,z 
2d p dlFj l 9 !\x,z)F^ pK (x,z 



-\d p dlp^(x,z)p$(x,z) 
-\d p dlp^(x,z)p%\x,z) 
-\d p dlp^{x,z)p^{x,z) 

-\d p dlpii\x,z)p( 9 ^(x,z) 

2d x ^ 9l ^(x,z)d z F^^(x,z) - 

2d p Flf{x,z)dlF^{x,z) - l -d p J${x,z)dlp^{x lZ ) 

2d p F^{x,z)diFW{x,z) - \d£pW^{x,z)dlpW{x,z) 

2d p F( 9l > K (x,z)diF^(x,z) - ±dZpW»«(x,z)dlpW(x,z) 

2d*F (9l)p Z{x,z)d z ! F^ p(3 {x,z) - ^d*pW p Z(x,z)d z pp( 9l ^(x,z) 

2d p F( 9l ^(x,z)d^ z F^f(x,z) - ldepW«{x,z)dj!pW{x,z)\ 

2d*F i9l)pl3 {x,z)d z fi F Wp Z{x,z) - ^p Wp(3 (x,z)dlp (9l)p ^x,z) 

2d p F(fJ\x,z)dlF^ pK (x,z) - ±d>ZpW(x,z)dlpW pK (x,zfj 

2^M^(x,z)S^W rf (s,z) - ^p^ 9l ^(x,z)d z ^ 9l ^(x,z) 

(A.l) 
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g 2 N 



a U a U (- 8F(f)aii ( x > z)p (f) ^(z, x) + 8p^ ap (x, z)fW«P(z, a)) 
+123£F« (x, z)df p (s) (x, z) + 12d p p^ (x, z)df F« (x, z) 
-l2d p x dlF {s \x,z)p i - s \x,z) - l2d p dlp ( - s \x,z)F^(x,z) 
+28lF^ h ) ( x , z)dPpW (z, x) - 2df />M ( Xj z )dPF^ h ) ( z , x ) 
-16 (d*d z pFWrf{x, z)p {9l) ^(x, z) + d*d z p P W p Z{x, z)FW)i#{x, z 



+8 
+8 
-4 
+8 
+8 
-4 
-4 
-4 
+16 



+4 
-8 
+4 
-8 
+4 
+4 



a/*d* F isi)pP( Xj z )pW^(x, z ) + d*d z p P w p P{x, z )f(*(x, 

d x d z F (9lW( Xj z )p i9l ^{x, Z) + d*d Z pp( 9l ^{x, z)FW(*( X , z)) 

9^|F^)^(x, *)pg>Or, z) + ^^(x, Z )^) (li z) ) 
d£dgFW"t( x , z)p f) (X) z) + dPdtpW"^, z)FW (x, z)) 
d£d%F$\x, z)pW»t( x , z) + pg } (x, z)FW»t( x , z)) 
dP%FW>{x, z)pW"«(x, z) + <9£dip«(x, z)F^)^(x, z)) 
d^diF^ix, z)pWP«( x , z) + <9£<9fp«(x, z)F^)^(x, z)) 
(d*FW p P{x, z)d z p p i9l) ^{x, z) + dfrpWWfa z)^F ( *(x, z)) 
^fS (x, z)df>pW>>t(x, z) + ^(x, z)d^F^(x, z)) 
d£FWf«{x, z)dlpW(x, z) + d£pWP«(x, z)9lFW(x, z)) 
d p F^(x, z)dlpl 9 ?(x, z) + dZpW>>«(x, z)%FW(x, z)) 
dPFWrt( x , z)d%pWtf(x, z) + d£pto l )(*(x, z)d z p F^^{x, z)) 
fr x FW»t{x, z)d?p%\x, z) + ^p^(x, z)dlFlf{x, z)) 
9JfW^(j, z)d%pW(<(x, z) + d*p( 9l ^{x, z)d%FW)Pt( x , z)) 
^FW(x, z)c\V^(x, z) + ^pg)(x, z)<9fF^(x, z)) 
9JF( j ^(x, z)d z ^ 9l ^(x, z) + d£pW>i*(x, z)^F« rf (i, z))" . 

(A.2) 



Fermion propagator corrections are shown in Figure 5 and yield 
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Figure 4: Corrections to gluon propagator. 





Figure 5: Corrections to fermion propagator. 



-2g 2 N 



-6 ( fW(z,x)F{x,z) + ]p[ f ){z,x)p(x,z) 



4 r rA 



(A.3) 



-2g 2 N 



(p^Or, Z )f(x, z) + z)p<$>(x, z)) 



+6 (p^ (z, x)F(x, z) - F ( J) (z, x)p(x, z)) 
Finally, the single ghost propagator correction is given in Figure 6, with 




(A.4) 



Figure 6: Correction to ghost propagator. 



-g 2 N 



2^fM( IiZ )^fW(x,2) - ±d£pW(x,z)%pW(x,z) 
+2FW(x,z)d£%FW( x , z ) - ± p W(x,z)dZ%p<* h )(x,z) 



(A.5) 



4 9h \x,z) = -2g 2 N[d^FW(x,z)d^ z p^ h \x,z) + d^pffd^F^(x,z) 

+F$\x, z)d%d»fi h \x, z) + p$\x, z)d£%FW( x , z)] . (A.6) 
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